WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 




PCT 

INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International Patent Classification 7 
H04J 3/06 



Al 



(11) International Publication Number: WO 00/16509 

(43) International Publication Date: 23 March 2000 (23.03.00) 



(21) International Application Number: PCT/US99/ 17848 

(22) International Filing Date: 5 August 1999 (05.08.99) 



(30) Priority Data: 

09/153,280 



15 September 1998 (15.09.98) US 



(71)(72) Applicant and Inventor: RAVIKANTH, Rayadurgam 
[IN/US]; 28-B, 26 Beacon Street, Burlington, MA 01803 
(US). 

(74) Agents: LYNCH, David, W. et al.; Altera Law Group, LLC, 
10749 Bren Road East, Opus 2, Minneapolis, MN 55343 
(US). 



(81) Designated States: AE, AL, AM, AT, AT (Utility model), AU, 
AZ, BA, BB, BG, BR, BY, CA, CH, CN, CR, CU, CZ, CZ 
(Utility model), DE, DE (Utility model), DK, DK (Utility 
model), EE, EE (Utility model), ES, FI, Fl (Utility model), 
GB, GD, GE, GH, GM, HR, HU, ID, 1L, IN, IS. JP, KE, 
KG, KP, KR, KZ. LC, LK, LR, LS, LT, LU, LV, MD, MG, 
MK, MN, MW, MX, NO, NZ, PL, PT, RO, RU, SD, SE, 
SG, SI, SK, SK (Utility model), SL, TJ, TM, TR, TT, UA, 
UG. UZ, VN, YU, ZA, ZW, ARIPO patent (GH, GM, KE, 
LS, MW, SD, SL, SZ, UG, ZW), Eurasian patent (AM, AZ, 
BY, KG, KZ, MD, RU, TJ, TM), European patent (AT, BE, 
CH, CY, DE, DK, ES, FI, FR, GB, GR, IE, IT, LU, MC, 
NL, PT, SE), OAPI patent (BF, BJ, CF, CG, CI, CM, GA, 
GN, GW, ML, MR, NE, SN, TD, TG). 



Published 

With international search report. 

Before the expiration of the time limit for amending the 
claims and to be republished in the event of the receipt of 
amendments. 



(54) Title: METHOD FOR ESTIMATING RELATIVE SKEW BETWEEN CLOCKS IN PACKET NETWORKS 



(57) Abstract 

A method, a system, and an arti- 
cle of manufacture for measuring a rel- 
ative skew between clocks in a packet 
based communication network. The 
method utilizes the timestamps avail- 
able in packets sent from a sender to 
a receiver, as well as the times at 
which the receiver receives such pack- 
ets. Such timestamps can be incorpo- 
rated in packets, or in a separate proto- 
col. 



200 





Get packet, N+1 












Get sender timestamp 
TN+i,sent,s, record receive 
timeTN+f.fece/vertD 







202 



204 



Compute delay: DelN* r.D 



206 



Compute Jitter 
J/v +t ,0 =De1/V+ i,D -De\N.D 



208 



J- 



Compute average ]Hter and average inter-packet time: 

EavtJN+T.DH2D+... .4jN+f,D)/N 
EaviJfU 1,sent t S -Tw,s«nf,S)=((T2,s*/if.S -T f ,sen/,s)+... .+(T/V+ 1,sent,S *TN t sent,S))fH 



210 



Compute relative skew: 
R/o/,5,D=Ea v(JW* f.0)/Ea v(TA/+ l.sent.S -T/V,*>nf,S)+1 



FOR THE PURPOSES OF INFORMATION ONLY 
Codes used to identify States party to the PCT on the front pages of pamphlets publishing international applications under the PCT. 



AL 


Albania 


ES 


Spain 


LS 


Lesotho 


SI 


Slovenia 


AM 


Armenia 


Ft 


Finland 


LT 


Lithuania 


SK 


Slovakia 


AT 


Austria 


FR 


France 


LU 


Luxembourg 


SN 


Senegal 


AU 


Australia 


GA 


Gabon 


LV 


Latvia 


sz 


Swaziland 


AZ 


Azerbaijan 


GB 


United Kingdom 


MC 


Monaco 


TD 


Chad 


BA 


Bosnia and Herzegovina 


GE 


Georgia 


MD 


Republic of Moldova 


TG 


Togo 


BB 


Barbados 


GH 


Ghana 


MG 


Madagascar 


. TJ 


Tajikistan 


BE 


Belgium 


GN 


Guinea 


MK 


The former Yugoslav 


TM 


Turkmenistan 


BF 


Burkina Faso 


GR 


Greece 




Republic of Macedonia 


TR 


Turkey 


BG 


Bulgaria 


HU 


Hungary 


ML 


Mali 


TT 


Trinidad and Tobago 


BJ 


Benin 


IE 


Ireland 


MN 


Mongolia 


UA 


Ukraine 


BR 


Brazil 


IL 


Israel 


MR 


Mauritania 


UG 


Uganda 


BY 


Belarus 


IS 


Iceland 


MW 


Malawi 


US 


United States of America 


CA 


Canada 


IT 


Italy 


MX 


Mexico 


uz 


Uzbekistan 


CF 


Central African Republic 


JP 


Japan 


NE 


Niger 


VN 


Viet Nam 


CG 


Congo 


KE 


Kenya 


NL 


Netherlands 


YU 


Yugoslavia 


CH 


Switzerland 


KG 


Kyrgyzstan 


NO 


Norway 


zw 


Zimbabwe 


ci 


Cflte d*I voire 


KP 


Democratic People's 


NZ 


New Zealand 






CM 


Cameroon 




Republic of Korea 


PL 


Poland 






CN 


China 


KR 


Republic of Korea 


PT 


Portugal 






CU 


Cuba 


KZ 


Kazakstan 


RO 


Romania 






CZ 


Czech Republic 


LC 


Saint Lucia 


RU 


Russian Federation 






DE 


Germany 


LI 


Liechtenstein 


SD 


Sudan 






DK 


Denmark 


LK 


Sri Lanka 


SE 


Sweden 






EE 


Estonia 


LR 


Liberia 


SG 


Singapore 







10 



PCT/US99/17848 

WO 00/16509 

METHOD FOR ESTIMATING RELATIVE SKEW BETWEEN CLOCKS IN 
Mtltiuurv, PACKET NETWORKS 

g APif fiROITNP OF THE INVENTION 

1 Field of the Invention . 

5 The present invention relates in general to communication systems, and more 

particularly, to a method and system for measuring/estimating the relative skew 
between clocks in a packet based communication networks. 

2 Description of Related Art. 

Estimating the relative skew between clocks on different host machines (or 
referred to network nodes) in a packet network environment, such as the Internet, is 
relatively new since real-time applications and usage based pricing are relatively 
new concepts. The issue of measuring one way transit times of packets is becoming 
increasingly important with the emergence of real-time applications such as voice 
and video over the Internet. This concept of one-way transit time measurements 

IS involves comparing (spe^^ 

todifferentclocks(e.g.thesenderandreceiverclocks). Two of the main sources of 

clocktimingdiscrepancyare^and^. Offsets can be readily corrected. The 
most widely used remedy for clock offsets is a form of synchronization (e.g. usmg 
Network Time Protocol (NTP), or Global Positioning System (GPS), etc.). It should 
20 benoted,*ough,*issynch^^ 

at one instant, but the clock skew will inherently cause the clocks to drift apart as 
time progresses. Accordingly, a clock skew can impose a serious clock timing 
discrepancy problem for a real-time application. 

Accordingly, there is a need to develop a method to detect/correct clock 
25 skews. Specifically.mereisaneedtodevelopamemodofmeasuring/estimatingthe 
relative skew between clocks on different host machines (or referred to network 
nodes) in a packet network environment such as the Internet. The present invention 
provides a solution to the above and other problems and offers advantages over the 
existing network systems. 
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SUMMARY OF THE INVENTION 
To overcome the limitations in the prior art described above, and to 
overcome other limitations that will become apparent upon reading and 
understanding the present specification, the present invention discloses a method and 
5 system for estimating the relative skew between clocks in a packet communication 
network. 

The present invention provides a method of detecting the existence of a clock 
skew between two different clocks. The clock offsets are predefined by a separate 
synchronization mechanism based on Network Time Protocol (NTP) and/or Global 
10 Positioning System (GPS) based systems. It is appreciated that other clock 

synchronization mechanisms can be used within the scope of the present invention. 

In one embodiment of the present invention, a method for determining a 
clock skew between at least a first and second clocks in a network includes: 
transmitting a plurality of packets from a first network node to a second network 
15 node, each of the packets having a timestamp which records time measured by the 
first clock at the first network node; receiving one of the packets at the second 
network node and recording time measured by the second clock at the second 
network node; obtaining the timestamp of the one packet recorded at the first clock; 
computing a time delay between a time of transmitting the one packet and a time of 
20 receiving the one packet; computing a jitter between the time delay for the one 
packet and the time delay for a next packet consecutively transmitted after the one 
packet; averaging jitters after transmitting the plurality of packets; averaging inter- 
packet times, each of the inter-packet times being a time interval between 
transmitting two consecutive of the plurality of packets at the first network node; 
25 and determining the clock skew by dividing an average of the jitters by an average of 

the inter-packet times. 

Further in one embodiment, the network is a packet based communication 

network. 

Still in one embodiment, the time recordings by the first and second clocks 

30 are relative to a global clock. 

The present invention is also a system of determining a clock skew between 
at least a first and second clocks in a network. In one embodiment, the system 
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includes means for transmitting a plurality of packets from a first network node to a 
second network node, each of the packets having a timestamp which records time 
measured by the first clock at the first network node; means for receiving one of the 
packets at the second network node and recording time measured by the second 
5 clock at the second network node; and means for determining a clock skew between 
the first and second clocks based on the timestamps and the time measured by the 
second clock. 

Further in one embodiment, the clock skew determining means includes means 
for obtaining the timestamp of the one packet recorded at the first clock; means for 
10 computing a time delay between a time of transmitting the one packet and a time of 
receiving the one packet; means for computing a jitter between the time delay for the 
one packet and the time delay for a next packet consecutively transmitted after the 
one packet; a first averager averaging jitters after transmitting the plurality of 
packets; and a second averager averaging inter-packet times, each of the inter-packet 
15 times being a time interval between transmitting two consecutive of the plurality of 
packets at the first network node; wherein the clock skew is determined by dividing 
an average of the jitters by an average of the inter-packet times. 

Further, the present invention provides an article of manufacture for a 
computer-based data processing system, the article of manufacture including a 
20 computer readable medium having instructions for causing a computer to perform a 
method for determining a clock skew between at least a first and second clocks in a 
network, includes: transmitting a plurality of packets from a first network node to a 
second network node, each of the packets having a timestamp which records time 
measured by the first clock at the first network node; receiving one of the packets at 
25 the second network node and recording time measured by the second clock at the 
second network node; obtaining the timestamp of the one packet recorded at the first 
clock; computing a time delay between a time of transmitting the one packet and a 
time of receiving the one packet; computing a jitter between the time delay for the 
one packet and the time delay for a next packet consecutively transmitted after the 
30 one packet; averaging jitters after transmitting the plurality of packets; averaging 
inter-packet times, each of the inter-packet times being a time interval between 
transmitting two consecutive of the plurality of packets at the first network node; 

3 
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and determining the clock skew by dividing an average of the jitters by an average of 

the inter-packet times. 

These and various other advantages and features of novelty which characterize 
the invention are pointed out with particularity in the claims annexed hereto and form a 
5 part hereof. However, for a better understanding of the invention, its advantages, and 
the objects obtained by its use, reference should be made to the drawings which form a 
further part hereof, and to accompanying descriptive matter, in which there are 
illustrated and described specific examples of an apparatus in accordance with the 
invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Referring now to the drawings in which like reference numbers represent 
corresponding parts throughout: 

Figure 1 illustrates a schematic view of one skew measurement scheme in a 
packet based communication network in accordance with the present invention. 
15 Figure 2 is a flow chart diagram of one skew measurement scheme to be 

implemented at a receiver or a destination network node in accordance with the 
present invention. 

r>FT AILED DESCRIPTION OF THE INVENTION 
20 The present invention relates in general to communication systems, and more 

particularly, to a method and system for measuring/estimating the relative skew 
between clocks in a packet based communication network. 

In the following description of the exemplary embodiment, reference is made 
to the accompanying drawings which form a part hereof, and in which is shown by 
25 way of illustration the specific embodiment in which the invention may be practiced. 
It is to be understood that other embodiments may be utilized as structural changes 
may be made without departing from the scope of the present invention. It is 
appreciated that the present invention can be used for real-time applications on 
Internet or other communication networks. 
30 The acronym used in the present invention are as follows: 

NTP Network Time Protocol 

GPS Global Positioning System 
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Time as measured by a global clock, which is considered to be 




thp 4t tnie time" 


Ti 


Timp n<: measured bv clock i. 


R 


ql-pw nf clock i with respect to the global clock 


5 Q 


offset of clock i with respect to the global clock at t g = 0. 
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QpnHpr or Source Network Node 


D 


Receiver or Destination Network Node 


Rs 


rinrV qVpw at S with respect to the global clock 


C s 


Clock offset at S with respect to the global clock at t g = 0. 


10 R D 


Clock skew at D with respect to the global clock 


C D 


Clock offset at D with respect to the global clock at t g = 0. 


Tfit.sent.S 


Time at which packet N is sent by a host machine of a sender 


S 

T Ni received,D Time at which packet N is received hy a host machine of a 


15 


receiver D 


Del KD 


Delay of packet N measured by the receiver D 


Del Ntt 


Delay measured by the global clock when R D =R S =1 and 



Jn+1,d Jitter measured by the receiver D 

20 Rrel,S,D Relative skew defined by Rd/Rs 

EMn+1,d> Average of the jitter measured by the receiver D. 

IPPM Internet Protocol Performance Metrics (IPPM is a working 



group 



of the Internet Engineering Task Force) 
25 rtp Real-time Transport Protocol 

The present invention provides a way of measuring/estimating the relative 
skew between different clocks in a packet based communication network. A 
preferred method utilizes timestamps available in packets sent from a sender S to a 
receiver D, as well as the times at which the receiver D receives such packets. The 
30 timestamps can be carried in the packet sent from the sender S to the receiver D, 
such that the timestamps can be read by the receiver's host computer. It is 

appreciated that the timestamps may be directly implemented in packets or 
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implemented by using a separate protocol so that the same measurement can be 
performed. Accordingly, a clock skew measurement scheme of the present 
invention can be integrated into a real-time application or can be used in conjunction 
with new protocols specifically designed to measure the performance of packet 
5 networks such as Internet 

Clock Issues 

As defined above, t g represents the time measured by a global clock, which 
is considered as the "true" time; Tj represents the time measured by a clock / at a 
10 network node i. The relation between 7/ and t g is as follows: 

Ti=Ri*t g +Q 

where R t is the skew of clock i with respect to the global clock, and Q is the offset 
of the clock i with respect to the global clock at tg=Q. In a preferred embodiment, 
the skew Rj is a constant parameter and does not vary with time. It is appreciated 
15 that the skew may vary with time in which cases they may be different. 

In general, clocks on different machines have different skews and offsets 
with respect to the global clock. As a result, there may be a difference between 
packet delays measured by the receiver D and the actual delay measured by the 
global clock. For clear illustration, a time point is defined in a network as follows: 
20 T N>sent ,s represents the time at which packet N was sent by the host machine of the 
sender S, and the time is measured by the clock of the sender S. Alternately, the 
second suffix can be received, and the third suffix can be any host index, such as g 
representing the global clock, or D representing the receiver's clock. 

From the above definition of times, 7jv,sen/,S =Rs* T N,sent,g +c S- 
25 Specifically, this is the timestamp (as, for instance, in the RTP timestamp) marked 
by S when the sender sends the packet N. Similarly, when the packet N is received 
at the destination node (i.e. the receiver D), the time that the packet N is received is 
T N ,received,D, (T N ,received,D = ^.received, g+Cn). Accordingly, the delay 
experienced by packet N, measured by the destination node D, is: 
30 DeI Ni o - Tn,received,D m TN,sent,s- 

It is appreciated that this measured delay coincides with the actual delay, 
Del Nig measured by the global clock, if and only if (i) Rd=Rs=\ and (ii) Cs=Cd. 
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The clock timing discrepancy are clock offset and skew. In the preferred 
embodiment, the present invention determines how much a relative offset between 
two clocks is. As mentioned in the background of the invention, offsets can be 
readily corrected by clock synchronization (using NTP, GPS etc.). It should be 
5 noted, though, that synchronization merely zeros the clock difference between two 
clocks at one moment, but the skew causes the clocks to drift apart as time 
progresses. Thus, the clock skew detection becomes an important issue in resolving 
clock timing discrepancy problem. A preferred method of detecting clock skew in 
accordance with the present invention is discussed below in details. 
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Clock Skew Detection 
One Of the main ideas behind a preferred method of detecting clock skews is 
that a sequence of packets transmitted from the source node (sender S) to the 
destination node (receiver D) are monitored. The delays of the packets are 
1 5 determined by the receiver D. It is appreciated that this paradigm can be used for 
one which uses a separate protocol to generate this sequence of packets, or for one 
which just uses data packets (time stamped UDP (User Datagram Protocol) packets, 
such as ones used by multimedia applications) to dynamically perform this 
measurement. 

20 As an example, consider the two consecutive packets N and N+l which are 

sent by the sender S to the receiver D. The delay suffered by the two packets are: 
Del N>D =T Nt received,D- T Ksen t iS mdDel N +l tD =T N +l t received,D- T N +l, S ent,S, 
respectively. Thus, a jitter measured by the receiver D is: 

Jn+1,D - Del N +j,o- Del N,D 
25 Substituting the detailed definitions of the quantities on the right hand side of 

this equality, J N +1,D = T N +i >re ceived,D- T N+l,sent,S- (JN,received,D- T Nl sent,s) 
=R D * T N +l,received,g + C D - ( R S * T N +i iS ent,g + - 

((R D * T N ,received,g + C D ) ■ ( R S * T N ,sent,g + Cs» 
=R D * T N +i ire ceived,g - ( R S * T N +l,sent,g) - 
30 «*D * T Nire ceived,g) - ( R S * T NtSen t,g)) 

Accordingly, the constant offset disappears when jitter is computed. To 
determine the relative skew Rrel s>£h the Rrel, s ,D = R c/ R S- 
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Jn+1,D = r D* (TN+l,received,g - TN+l,sent,g • <JKreceived,g - TN,sent,g)) 
+ (R D -Rs> * (T H+l,sent,g - T N>S ent,g) 
= R D * (Del N + i ng - Del Nig ) + (R D - RJ * (TN+l,sent,g - TN.sent.g) 
= R D * (JN+J,g) + (Rd-As) * (TN+l,sent,g - TN.sent.g). 
5 where J N + j ig is the jitter measured by the global clock. T N +i iSe „t,g and T NtSen t t g 
can then be rewritten in terms of the timestamp put in by the senders. Thus, 

Jn+1,d = Rd * (JN+l,g) + (*D-*st * ((TN+l,sent,S- Cs)/R s - <T N>S ent,s- 
Cs)/Rs) 

= RD(JN+l.g) + (TN+lsent.S- T N , S era,S> * (Rd- r s)^S 
10 =RDVN+l,g) + (TN+l.sent,S- T N ,sent,S> * (Rrel s>D -V- 

Accordingly, the above expression allows one to measure the relative skew 
Rrel s d- Q y taking averages (ensemble or time averages, either denoted by EavQ) 
on both sides of the above equality, then 

E av (J N +l,D) =Eav(T N +l,sent,S • T KS ent,s) * (Rrel, SiDr l), 
1 5 where Eav(J N +i ig ) = 0. The validity of this fact can be seen in the stochastic sense 
since a jitter is the difference of two delays, which may be assumed to be identically 
distributed in the stochastic sense. In the time average sense, the explanation for this 
is that a non-zero mean can imply that the delays are exploding to infinity (positive 
average case) or they are going to be zero (negative average case). Now, note that if 
20 the source node S is sending packets at some average time interval denoted by 
Tinter-packeh th en 

EovVn+I.d) = Ti„ter-packet*(R™k,D-V- 
In the above relation, all quantities, except for the relative skew Rrel s ,D, are 
known. Therefore, the relative skew can be readily dividing the average jitter 
25 Eav(J N +l, D ) by the average time interval Ti nte r-packet- In a preferred embodiment, 
the time average approach is used. It is appreciated that other average approach, 
such as ensemble average approach, can be used within the scope of the present 
invention. 

30 Usage of the skew determination test 

In accordance with the above scheme, in one embodiment, a sequence of 
packets are sent from a sender to a receiver. A number of observations are made to 
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obtain a more accurate estimate of the clock skew. The packets can be sent at a 
fixed intervals of time. It is appreciated, however, that the packets need not be sent 
at fixed intervals of time since all it is needed is the time average of the timestamps 
from the sender. Thus, if, for example, a Poisson sampling is used to generate the 

5 packet generation times, this method can remain unchanged. 

The above method can be utilized in conjunction with existing applications. 
Real-time applications, such as some video-conferencing applications, e.g. H.323 (a 
protocol which enables packet based multimedia communication system) based 
video-conferencing applications, may need to know timing information, and hence 

1 0 may require the receiver to estimate the relative skew between its clock and the 
sender's clock. In such a situation, the timestamps available in the Real-time 
Transport Protocol (RTP) header of the data packets may be sufficient for the 
receiver to estimate/measure the relative skew. Thus, this method can potentially be 
integrated into applications. 

15 

Relative skew measurement implementation 
One implementation of the above clock skew determination method is 
described below. 

Figure 1 illustrates a block diagram of skew measurement scheme in a packet 
20 based communication network 100. Packets 102, each having a timestamp 104 
(either incorporated in an application or in a separate protocol), are sent from a 
source node 106 to a destination node 108 over the network 100. The destination 
node D 108 measures the relative skew between a clock 1 10 at the source node 106 
and a clock 1 12 at the destination node 108 by using a stream of consecutive packets 
25 from the source node 106 to the destination node 108. 

The determination of the clock skew is made at the destination node 108 
through the use of the relations specified in the equation described before. Figure 2 
shows an operation of the relative skew measurement at the destination node 108. 
The computation of a clock skew is made at the reception of each packet. The 
30 estimate of the overall clock skew becomes more and more accurate as more and 
more packets are received. It is appreciated that the flow chart shown in Figure 2 is 
only for illustration, not for limitation. For example, the sequence of the 
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measuring/estimating steps can be altered within the principles of the present 
invention. 

In Figure 2, a packet N+l is received by a receiver in box 200. The receiver 
obtains the sender timestamp from the packet meanwhile recording the receiving 

5 time in box 202. A delay Del N+UD is determined in box 204. A jitter is then 

computed in box 206. In box 208, an average jitter and an average inter-packet time 
are computed by averagers which are commonly used by a person skilled in the art. 
Finally, a relative skew Rrel SD is determined in box 210. 

Accordingly, the present invention provides a method for 

1 0 measuring/estimating the relative clock skew between network nodes in a packet 
based communication network. The potential use of such a scheme can be ranged 
from integration into applications to use by network operators as a performance 
measurement technique. Further, it is appreciated that the above clock skew 
measurement technique can be generally used for many other clock skew correction 

15 applications. 

The foregoing description of the exemplary embodiment of the invention has 
been presented for the purposes of illustration and description. It is not intended to 
be exhaustive or to limit the invention to the precise form disclosed. Many 
modifications and variations are possible in light of the above teaching. It is 
20 intended that the scope of the invention be limited not with this detailed description, 
but rather by the claims appended hereto. 
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WHAT IS CLAIMED IS: 

1 . A method for determining a clock skew between at least a first and second 
clocks in a network, comprising: 

transmitting a plurality of packets from a first network node to a second 
5 network node, each of the packets having a timestamp which records time measured 
by the first clock at the first network node; 

receiving one of the packets at the second network node and recording time 
measured by the second clock at the second network node; 

obtaining the timestamp of the one packet; 
1 0 computing a time delay between a time of transmitting the one packet and a 

time of receiving the one packet; 

computing a jitter between the time delay for the one packet and the time 
delay for a next packet consecutively transmitted after the one packet; 

averaging jitters after transmitting the plurality of packets; 
15 averaging inter-packet times, each of the inter-packet times being a time 

interval between transmitting two consecutive of the plurality of packets at the first 
network node; and 

determining the clock skew by dividing an average of the jitters by an 
average of the inter-packet times. 

20 2. The method of claim 1 , wherein the network is a packet based 
communication network. 

3. The method of claim 1 , wherein time recordings by the first and second 
clocks are relative to a global clock. 

4. A system of determining a clock skew between at least a first and second 
25 clocks in a network, comprising: 

a transmitter for transmitting a plurality of packets from a first network node 
to a second network node, each of the packets having a timestamp which records 
time measured by the first clock at the first network node; 
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a receiver for receiving one of the packets at the second network node and 
recording time measured by the second clock at the second network node; and 

means for determining a clock skew between the first and second clocks 
based on the timestamps and the time measured by the second clock. 

5 5. The system of claim 4, wherein the clock skew determining means includes: 
means for obtaining the timestamp of the one packet recorded at the first 

clock; 

means for computing a time delay between a time of transmitting the one 
packet and a time of receiving the one packet; 
1 0 means for computing a jitter between the time delay for the one packet and 

the time delay for a next packet consecutively transmitted after the one packet; 

a first averager averaging jitters after transmitting the plurality of packets; 

and 

a second averager averaging inter-packet times, each of the inter-packet 
15 times being a time interval between transmitting two consecutive of the plurality of 
packets at the first network node; wherein the clock skew is determined by dividing 
an average of the jitters by an average of the inter-packet times. 

6. The system of claim 4, wherein the network is a packet based 
communication network. 

20 7. The system of claim 4, wherein the time recordings by the first and second 
clocks are relative to a global clock. 

8. An article of manufacture for a computer-based data processing system, the 
article of manufacture comprising a computer readable medium having instructions 
for causing a computer to perform a method for determining a clock skew between 
25 at least a first and second clocks in a network, comprising: 

transmitting a plurality of packets from a first network node to a second 
network node, each of the packets having a timestamp which records time measured 
by the first clock at the first network node; 
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receiving one of the packets at the second network node and recording time 
measured by the second clock at the second network node; 
obtaining the timestamp of the one packet; 

computing a time delay between a time of transmitting the one packet and a 
5 time of receiving the one packet; 

computing a jitter between the time delay for the one packet and the time 
delay for a next packet consecutively transmitted after the one packet; 

averaging jitters after transmitting the plurality of packets; 

averaging inter-packet times, each of the inter-packet times being a time 
1 0 interval between transmitting two consecutive of the plurality of packets at the first 
network node; and 

determining the clock skew by dividing an average of the jitters by an 
average of the inter-packet times. 



9. The article of claim 8, wherein the network is a packet based communication 
15 network. 

10. The article of manufacture of claim 8, wherein time recordings by the first 
and second clocks are relative to a global clock. 
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200 



Get packet ,N+1 



202 



Get sender timestamp 
TN+i,sent,S, record receive 
time TN+1,received,D 



Compute delay: DelA/+/,D 

i 



204 



206 



Compute Jitter: 
JN+1.D =DelN+ 1,D -Delw.D 



208 



Compute average jitter and average inter-packet time: 
Eav(JA/+f,D)=j2D+... .+JW+J,D)/N 
Eav{TN+1,sent,S-TN,sent,S)={{T2,sent,S-Tl,sent,s)+... .+{TN+1,sent l S-TN,sent,S))/N 



210 



Compute relative skew: 
Rre/,S.D=Ea v(jN+f ,0)/Eav(TA/+ 1,sent,S •TN.sent.sW 
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